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SUMMARY 


Ten years ago, computer simulations of the flow field around the 
Space Shuttle orbiter were limited to inviscid calculations for 
the windward side of the forebody and viscous calculations for 
selected two-dimensional problems. Advances in computer hardware 
and numerical methods during the past 10 years have made it pos- 
sible to calculate viscous flow over the complete orbiter con- 
figuration at angle of attack. The equations solved are the 
Reynolds-averaged, Navier-Stokes equations, simplified by either 
the thin-layer or parabolized approximation. An algebraic eddy 
viscosity model is used for turbulent flow. The free stream is 
assumed to be a perfect gas for wind tunnel conditions and a real 
gas in thermodynamic equilibrium for flight conditions. Four ex- 
amples of recent computer simulations will be presented. Flow 
field results include oil flow patterns on the surface and Mach 
number contours, isobars, and cross-flow velocity vectors in the 
shock layer. 


INTRODUCTION 


Preliminary computer si- 
mulations of the flow field 
around the Space Shuttle or- 
biter, shown in Fig. 1, were 
presented at the 12th Interna- 
tional Symposium on Space 
Technology and Science held in 
Tokyo in 1977 (Ref. 1). The 
results were limited to invis- 
cid calculations for the wind- 
ward side of the forebody and 
viscous calculations for two- 
dimensional problems. Advan- 
ces in computer hardware and 
numerical methods have made it 
possible to calculate viscous 
flow over the complete orbiter 
configuration at angle of at- 
tack. 



Fig. 1 Space Shuttle orbiter 


Ten years ago the Control Data Corporation 7600 with a peak 
processing speed less than 10 million floating point 
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operations/sec and a main memory of half a million words was the 
largest computer generally available. The current generation of 
supercomputers was just being introduced to the scientific com- 
munity, and NASA Ames Research Center had the unique ILLIAC IV. 
Today, supercomputers are available with peak speeds exceeding 
one billion floating point operations/sec and with main memories 
of 250 million words. Their availability extends to academia, 
industry, and governmental organizations worldwide as shown in 
Fig. 2. Advances have also been made in peripheral equipment for 
storing data and displaying results of computations. 



2 U.S. 

2 Japan 


Fig. 2 Worldwide Distribution of Supercomputers 


Concurrent with advances in computer speed and memory, 
faster numerical methods have been developed to drastically 
reduce the computer time required to solve complex flow fields. 
Ten years ago, explicit finite difference methods, which require 
small time steps to achieve stability, were used, and the 
MacCormack method was the most popular. Today, implicit meth- 
ods, which are not so restricted, are used in many computer 
codes, and the Beam-Warming class of methods has replaced the 
MacCormack method in popularity. Advances in grid generation 
techniques also make it possible to produce physically appropri- 
ate computational meshes around complex aircraft configurations. 

Four examples of recent computer simulations will be shown 
for viscous flow over the Space Shuttle orbiter. The equations 
solved are the Reynolds-averaged, Navier-Stokes equations, sim- 
plified by either the thin-layer or parabolized approximation to 
facilitate their solution. For the high Reynolds numbers nor- 
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mally associated with aerodynamic configurations, the viscous ef- 
fects are confined to the neighborhood of the body, which results 
in the derivatives of the viscous terms normal to the surface 
being much larger than the corresponding derivatives along the 
surface. When the latter are neglected, the thin-layer approxi- 
mation results. If the flow is also steady and supersonic, the 
parabolized Navier-Stokes equations can be obtained by simpli- 
fying the streamwise pressure gradient within the viscous layer. 
In this case the resulting equations may be solved by a space- 
marching method, proceeding downstream step by step from one 
cross-sectional plane to another. Available computer storage can 
be used effectively for this approach because just a fraction of 
the flow field properties need be stored at any given time. 
The solution algorithms used in the following examples are all 
based on the implicit Beam-Warming method. 


RESULTS 


Example 1 


The first example is a solution by Chaussee, Rizk, and 
Buning (Ref. 2) for the parabolized Navier-Stokes equations with 
a perfect gas free stream representative of wind tunnel condi- 



Y 

(a) Computational grid 


Fig. 3 Perfect gas solution, 
(Ref. 2) 



Y 

(b) Mach contours 


M„ = 7.9, a = 25°, x/L =0.66 
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tions. The free-stream Mach number, M„ , was 7.9, the angle of 
attack, a, was 25°, and the Reynolds number based on vehicle 
length was 1.37 million. Turbulent flow was assumed over the 
surface with the eddy viscosity calculated from the Baldwin-Lomax 
algebraic model (Ref. 3) . The calculations were initiated at the 
cross section located at x/L = 0.0522, where x is the axial 
distance measured from the nose and L is the total length of the 
orbiter. The axial flow was supersonic at this cross section, 
and flow properties were obtained from a blunt body code. The 
flow field properties were then marched downstream as far as the 
cross section, x/L = 0.66, where the interaction of the bow 
shock wave with the wing shock caused a local subsonic region 
which invalidated the solution procedure. The computational grid 
at x/L = 0.66 is shown in Fig. 3a. The outer boundary is the 
bow shock wave, and points are clustered near the surface to re- 
solve the viscous layer. Mach number contours are shown in 
Fig. 3b, where coalescence of contours indicates the locations of 
the wing shock and cross flow shocks on the wing and upper body. 
Simulated oil flow patterns over the upper orbiter surface are 
shown in Fig. 4. Separated flow is indicated on the strake-wing 
and payload bay door on the upper half of the fuselage. 



Fig. 4 Oil flow pattern, M m = 7.9, a = 25° (Ref. 2) 


Example 2 

The second example is a solution by Rizk, Steger , and 
Chaussee (Ref. 4) for the thin-layer Navier-Stokes equations with 
a perfect gas free stream. The flow field was subdivided into 
four regions, as shown in Fig. 5, to accommodate different grids 
and numerical methods and to make effective use of the available 
computer storage. Region I was solved first, and results at a 
downstream plane, where the flow was supersonic, was used as in- 
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put for Region II calcula- 
tions. This procedure was re- 
peated for Regions III and IV. 
Time-dependent procedures were 
used in all four regions, al- 
though a space-marching method 
could have been implemented in 
Region III to conserve com- 
puter time. 

Results are presented at 
the aft fuselage station, 
x/L = 0.91, in Fig. 6 for 
Moo = 1.4, a = 0°, and turbu- 
lent flow, using the Baldwin- 
Lomax model. Details of the 
computational grid near the 
body are shown in Fig. 6a, in- 
cluding the wing, orbital 
maneuvering system pod, and 
vertical tail. Isobars are 
shown in F ig . 6b . 



Fig. 5 Computational regions 
(Ref. 4) 
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(a) Computational grid 

Fig. 6 Perfect gas solution, M 
(Ref. 4) 

Example 3 

The third example is a solu 



(b) Isobars 

= 1.4, a = 0°, x/L = 0.91 


ion by Prabhu and Tannehill 
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(Ref. 5), which incorporates real-gas effects in a parabolized 
Navier-Stokes code. The thermodynamic properties were obtained 
from curve fits to tables that are applicable for equilibrium air 
at high temperatures, including the effects of dissociation and 
ionization. A modified orbiter geometry, deleting the canopy, 
vertical tail, and orbital maneuvering system pods, was used for 
the calculations. 

The parabolized Navier-Stokes equations were solved, begin- 
ning at x/L = 0.015, with initial conditions obtained from a 
blunt-body code for the nose region. The flow field properties 
were marched downstream until terminated at x/L = 0.40 by numeri- 
cal difficulties attributed to the grid generation method. Mach 
number contours at x/L = 0.40 are shown in Fig. 7 for laminar 
flow at M« = 21, a = 40°, and an altitude of 71.3 km. 

Example 4 

The last example is a solution by Balakrishnan (Ref. 6), 
which also incorporates real-gas effects in a parabolized Navier- 
Stokes code. In this case, however, the species concentrations 
and thermodynamic properties were calculated at each grid point 


400 



Y, in. 

Fig. 7 Real gas solution, Mach Fig. 8 Real gas solution, 
contours, M„, = 21, a = 40°, cross-flow velocity vectors, 

altitude = 71.3 km, x/L = 0.4 M„ = 13, a = 0°, alti- 

(Ref. 5) tude = 55.8 km, x/L = 0.53 

(Ref. 6) 
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using formulations obtained from the kinetic theory of gases. 
Air was assumed to be composed of seven species in thermodynamic 
equilbrium, including oxygen, nitrogen, and nitric oxide mole- 
cules; oxygen and nitrogen atoms; nitric oxide ions; and elec- 
trons. Although requiring more computer time, this approach 
provides species mole fractions necessary for calculating radia- 
tive heat transfer and is readily extendible to nonequilibrium 
flow calculations. The simplified geometry of the third example 
was also used for these calculations. 

Initial conditions at the cross section, x/L = 0.026, which were 
obtained from a blunt body solution, were marched downstream to 
x/L = 0.61. Cross-flow velocity vectors at the cross section, 
x/L = 0.53, are shown in Fig. 8 for laminar flow at = 13, 
a = 0°, and an altitude of 55.8 km. Flow reversal is indicated 
at the wing-body juncture. 


CONCLUDING REMARKS 


Computer simulations can now provide detailed information 
about the viscous flow around the space shuttle orbiter for a 
wide range of free-stream Mach and Reynolds numbers and gas 
properties. Computer codes validated by wind tunnel tests can be 
applied at flight conditions which cannot be reproduced in test 
facilities. The aerodynamic design of the Space Shuttle orbiter 
was fixed 10 years ago, but computer simulations are useful 
today, for example, for studies of a crew escape system during 
the launch phase. In addition, the computer codes can be used 
for preliminary design studies of advanced aerospace vehicles. 
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